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o Abstract 



In awake animals, the activity of the cerebral cortex is highly complex, with neurons firing iiTcgularly with 
apparent Poisson statistics. One way to characterize this complexity is to take advantage of the high inter- 
connectivity of cerebral cortex and use intracellular recordings of cortical neurons, which contain information 
about the activity of thousands of other cortical neurons. Identifying the membrane potential (V,„) to a stochas- 
tic process enables the extraction of important statistical signatures of this complex synaptic activity. Typically, 
one estimates the total synaptic conductances (excitatory and inhibitory) but this type of estimation requires at 
least two V,„ levels and therefore cannot be applied to single V,„ traces. We propose here a method to extract ex- 
citatory and inhibitory conductances (mean and variance) from single V,„ traces. This "VmT method" estimates 
conductance parameters using maximum likelihood criteria, under the assumption that synaptic conductances 
are described by Gaussian stochastic processes and are integrated by a passive leaky membrane. The method is 
illustrated using models and is tested on guinea-pig visual cortex neurons in vitro using dynamic-clamp exper- 
iments. The VmT method holds promises for extracting conductances from single-trial measurements, which 
has a high potential for in vivo apphcations. 



In awake animals, neurons in different cortical structures display highly irregular spontaneous firing (Evarts, 
1964; Steriade and McCarley, 1990). This level of firing must be considered together with the dense connec- 
tivity of cerebral cortex. Each pyramidal neuron receives between 5000 and 60000 synaptic contacts and a 
large part of this connectivity originates from the cortex itself (DeFelipe and Farinas, 1992; Braitenberg and 
SchiJz, 1998). As a consequence, many synaptic inputs are simultaneously active onto cortical neurons in intact 
networks. Indeed, intracellular recordings in awake animals reveal that cortical neurons are subjected to an 
intense synaptic bombardment and, as a result, are depolarized and have a low input resistance (Matsumara 
et al., 1988; Baranyi et al., 1993; Steriade et al., 2001) compared to brain slices kept in vitro. This activity 
is also responsible for a considerable amount of subthreshold fluctuations, called "synaptic noise". This noise 
level and its associated "high-conductance state" greatly affect the integrative properties of neurons (reviewed 
in Destexhe et al., 2003; Destexhe, 2007; Longtin, 2008). 

To characterize this activity, one must use appropriate methods to measure not only the mean conductance level 
of excitatory and inhibitory inputs, but also their level of fluctuations which is quantified by the conductance 
variance. Many methods exist to extract mean conductances (see review by Monier et al., 2008), but few meth- 
ods have been proposed to extract the variances. One class of methods consists of identifying the membrane 
potential (V,„) to a multidimensional stochastic process (Rudolph and Destexhe, 2003c, 2005). The so-called 
"VmD method" (Rudolph et al., 2004) enables the extraction of mean and variance of excitatory and inhibitory 
conductances. This method, however, requires to use at least two different levels of V,,, activity, which prevents 
its application to single-trial measurements. 

More recently, a related method was proposed to extract spike-triggered average conductances from the V,„ 
(Pospischil et al., 2007). In this case, one needs to extract synaptic conductances time courses from a single 
trace (the spike-triggered V^), which was made possible using a maximum likelihood method. 

In the present paper, we attempt to merge these concepts and provide a method to estimate the excitatory 
and inhibitory conductances and their variances. This "VmT" method provides similar estimates as the VmD 
method, but it does so by using a maximum likelihood estimation, and thus can be applied to single V„, traces. 



Experimental procedures 



Models 

All methods and analyses shown in this paper are based on a stochastic model of synaptic background activity, 
which was modeled by fluctuating conductances (Destexhe et al., 2001). In this model, the synaptic conduc- 
tances are stochastic processes, which in turn influence V„j dynamics. According to this "point-conductance" 
model, the membrane potential dynamics is described by the following set of equations: 
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where C denotes the membrane capacitance, Ig^t a stimulation cuiTcnt, Gi the leak conductance and Ei the leak 
reversal potential. ge{t) and gi{t) are stochastic excitatory and inhibitory conductances with respective reversal 



potentials Ee and £",. The excitatory synaptic conductance is described by Ornstein-Uhlenbeck stochastic pro- 
cesses (Eq. |2l), where g^o and d], are, respectively, the mean value and variance of the excitatory conductance, 
Tg is the excitatory time constant, and i,e{t) is a Gaussian white noise source with zero mean and unit standard 
deviation. The inhibitory conductance gi{t) is described by an equivalent equation (Eq.O with parameters gio, 
oj, Xi and noise source ^,(f)- For details about Ornstein-Uhlenbeck stochastic processes, see the review by 
Gillespie (1996). 

This model was simulated together with a leaky integrator mechanism with no spiking mechanism because 
the method only applies to subthreshold activity. All numerical simulations were performed in the NEURON 
simulation environment (Hines and Carnevale, 1997) and were run on PC -based workstations under the LINUX 
operating system. 



Experiments 

In vitro experiments were performed on 0.4 mm thick coronal or sagittal slices from the lateral portions of 
guinea-pig occipital cortex. Guinea-pigs, 4-12 weeks old (CPA, Olivet, France), were anesthetized with sodium 
pentobarbital (30 mg/kg). The slices were maintained in an interface style recording chamber at 33-35 °C. 
Slices were prepared on a DSK microslicer (Ted Fella Inc., Redding, CA) in a slice solution in which the NaCl 
was replaced with sucrose while maintaining an osmolarity of 307 mOsm. During recording, the slices were 
incubated in slice solution containing (in mM): NaCl, 124; KCl, 2.5; MgS04, 1.2; NaHP04, 1.25; CaCb, 2; 
NaHCOs, 26; dextrose, 10, and aerated with 95% O2, 5% CO2 to a final pH of 7.4. Intracellular recordings 
following two hours of recovery were performed in deep layers (layer IV, V and VI) in electrophysiologically 
identified regular spiking and intrinsically bursting cells. Electrodes for intracellular recordings were made on a 
Sutter Instruments P-87 micropipette puller from medium-walled glass (WPI, IBFIOO) and beveled on a Sutter 
Instruments beveler (BV-IOM). Micropipettes were filled with 1.2 to 2M potassium acetate and had resistances 
of 60-100 MD. after bevehng. 

All research procedures concerning the experimental animals and their care adhered to the American Phys- 
iological Society's Guiding Principles in the Care and Use of Animals, to the European Council Directive 
86/609/EEC and to European Treaties series no. 123, and was also approved by the local ethics committee 
"Ile-de-France Sud" (certificate no. 05-003). 

The dynamic-clamp technique (Robinson et al., 1993; Sharp et al., 1993) was used to inject computer-generated 
conductances in real neurons. Dynamic-clamp experiments were run using the hybrid RT-NEURON environ- 
ment (Sadoc et al., 2009), which is a DSP-based system using a modified version of the NEURON simulation 
environment (Hines and Carnevale, 1997). The dynamic-clamp protocol was used to insert the fluctuating 
conductances underlying synaptic noise in cortical neurons using the point-conductance model, similar to a 
previous study (Destexhe et al., 2001). According to Eq. ([T]) above, the injected cuiTcnt is determined from the 
fluctuating conductances ge{t) and gi{t) as well as from the difference of the membrane voltage from the respec- 
tive reversal potentials, loynciamp = —ge {V — Ee) — gi [V — Eij . The contamination of the measured membrane 
voltage V by electrode artifacts due to simultaneous cunent injection was avoided through the use of Active 
Electrode Compensation (AEC), a novel, high resolution digital on-line compensation technique (Brette et al. 
2008). 



Results 



We present here a method for extracting conductances by associating the V,„ to a stochastic process, and which 
uses maximum likelihood criteria to perform this extraction from single V^ traces. We first describe the method, 
then test it successively using computational models and real neurons in vitro. 



The VmT method: extracting conductances from single V„j traces 

The starting point of the method is to search for the "most likely" conductances parameters ( g^o, gm, Oe and 
Oi) that are compatible with an experimentally observed Y,„ trace. We start from the point-conductance model 
(Eqs.[T]-[3l), which is discretized in time with a step-size At. Eq.[T]can then be solved for ^^, which gives: 



C JV'--£'l g^.iV'-E,) V'^+^-V' 
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where Tl = C/Gl. 

Since the series V'^ for the voltage trace is known, g^ has become a function of g^. In the same way, we solve 
Eqs.|2]-[3]for ^J, which have become Gaussian distributed random numbers. 
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where s stands for e, /. 



There is a continuum of combinations {ge^^gf^^} that can advance the membrane potential from V^^^ to 
V^'^'^, each pair occurring with a probability 

p' := pi!A*\>sr\£J,) = ^e-¥.if*9''. = l-e-i!''', (6) 
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These expressions are identical to those derived previously for calculating spike-triggered averages (Pospischil 
et al., 2007), except that no implicit average is assumed here. 

Thus, to go one step further in time, a continuum of pairs (gg^\g,-^^) is possible in order to reach the (known) 
voltage V^+^. The quantity p*^ assigns to all such pairs a probability of occurrence, depending on the previous 
pair, and the voltage history. Ultimately, it is the probability of occurrence of the appropriate random numbers 
^^' and l^^j that relate the respective conductances at subsequent time steps. It is then straightforward to write 
down the probability p for certain conductance series to occur, that reproduce the voltage time course. This is 
just the probability for successive conductance steps to occur, namely the product of the probabilities p'^: 

P = Y\p\ (8) 

k=Q 

given initial conductances g^, g°. However, again, there is a continuum of conductance series 

{ge^Si}l=\,...,n+\, 



that are all compatible with the observed voltage trace. We define a likelihood function /( V*^ , 6) , 6 = (g^o , gio , ^e , (^i 
that takes into account all of them with appropriate weight. We thus integrate Eq.[8]over the unconstrained con- 
ductances gg and normalize by the volume of configuration space: 



m=UsUsip' 
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where only in the nominator g^ has been replaced by Eq.|4] This expression reflects the likelihood that a specific 
voltage series {V*"} occurs, normalized by the probability, that any trace occurs. The most likely parameters 
6 giving rise to {V'^} are obtained by maximizing (or minimizing the negative of) /(V*^,0) using standard 
optimization schemes (Press et al., 2007). 



Test of the VmT method using model data 

We tested the method in detail in its applicability to voltage traces that were created using the same model (leaky 
integrator model). To this end, we performed simulations scanning the (gt.o,g;o)-plane, and subsequently tried 
to re-estimate the conductance parameters used, but only from the V„, activity. For each parameter set {geOigio, 
Oe, <^i) the method was applied to ten samples of 5000 data points (corresponding to 250 ms each) and the 
average was taken subsequently. The conductance standard deviations (SDs) were chosen to be one third of 
the respective mean values, other parameters were assumed to be known during re-estimation (C = 0.4 nF, 
gL = 13.44 nS, El = -80 mV, E^ = mV, E, = -75 mV, x^ = 2.728 ms, X; = 10.49 ms), the time step was 
dt = 0.05 ms. Also, we assumed that the total conductance gtot was known - this is the inverse of the apparent 
input resistance, which is generally known. This assumption is not mandatory, but the estimation becomes 
more stable. The likelihood function given by Eq. |9]was thus only maximized with respect to geo,<^e and a,. 
Fig. [T] summarizes the results obtained. 

The mean conductances are well reproduced over the entire scan region. An exception is the estimation of g,o in 
the case where the mean excitation exceeds inhibition by several-fold, a situation which is rarely found in real 
neurons. The situation for the SDs is different. While the excitatory SD is reproduced very well in the whole 
ai^ea under consideration, this is not necessarily the case for inhibition. Here, the estimation is good for most 
parts of the scanned region, but shows a considerable deviation along the left and lower boundaries. These are 
regions where the transmembrane cuiTcnt due to inhibition is weak, either because the inhibitory conductance 
is weak (lower boundary) or because it is strong and excitation is weak (left boundary), such that the mean 
voltage is close to the inhibitory reversal potential and the driving force is small. In these conditions, it seems 
that the effect of inhibition on the membrane voltage cannot be distinguished from that of the leak conductance. 

This point is illustrated in Fig. |2] The relative deviation between a, in the model and its re-estimation depends 
on the ratio of the transmembrane cunent due to inhibitory (/,) and leak (7^) conductance. The estimation 
fails when the inhibitory current is smaller or compai^able to the leak cun^ent, but it becomes very reliable 
as soon as the ratio /,///, becomes larger than 1.5-2. Some points, however, have large en^ors although with 
dominant inhibition. These points have strong excitatory conductances (see gray scale) and correspond to the 
upper right corner of Fig.[ljD. The error is due to aberrant estimates for which the predicted variance is zero; 
in principle such estimates could be detected and discarded, but no such detection was attempted here. Besides 
these particular combinations, the majority of parameters with strong inhibitory conductances gave acceptable 
errors. We conclude that the estimation of the conductance variances will be most accurate in high-conductance 
states, where inhibitory conductances are strong and larger than the leak conductance. 



Effect of recording noise 

The unavoidable presence of recording noise may present a problem in the application of the method to record- 
ings from real neurons. Fig. [3] (left) shows how low-amplitude white noise added to the voltage trace of a leaky 
integrator model impairs the reliability of the method. A Gaussian-distributed white noise was added to the 
voltage trace at every time step, scaled by the amplitude given in the abscissa. Different curves con^espond to 
different pairs (geo, gio) colored as a function of the total conductance. The noise has an opposite effect on the 
estimation of the conductance mean values. While the estimate of excitation exceeds the real parameter value, 
for inhibition the situation is inverted. However, one has to keep in mind that both parameters are not estimated 
independently, but their sum is kept fixed. In contrast, the estimates for the conductance SDs always exceed 
the real values, and they can deviate by almost 500% for a noise amplitude of 10 juY. Here, the largest errors 
generally correspond to lowest conductance states. Clearly, in order to apply the method to recordings from 
real neurons, one needs to restrain this noise sensitivity. 

Fortunately, this noise sensitivity can be reduced by standai^d noise reduction techniques. For example, pre- 
processing and smoothing the data using a Gaussian filter greatly diminishes the amplitude of the noise, and 
consequently improves the estimates according to the new noise amplitude (see Fig.[3l right panels). Too much 
smoothing, however, may result in altering the signal itself, and may introduce enors. It is therefore preferable 
to use smoothing at very short time scales (SD of 1-4 data points, depending on the sampling rate). In the next 
sections, we preprocessed the experimental voltage traces with a Gaussian filter with a SD of 3 data points. 



Application of the VmT method to in vitro data 

We next tested the method on in vitro recordings (n=5 cortical neurons, 17 injections) using dynamic clamp 
experiments (Fig. |4]). As in the model, the stimulus consisted of two channels of fluctuating conductances 
representing excitation and inhibition. The conductance injection spanned values from low-conductance (of the 
order of 5-10 nS) to high-conductance states (50-160 nS). It is apparent from Fig.|4]that the mean values of the 
conductances (geo, gio) ai"e well estimated, as expected because the total conductance is known in this case. 

However, the estimation is subject to larger enors for the standard deviations of the conductances (a^, a,)- 
In addition, the eiTor on estimating variances is also linked to the accuracy of the estimates of synaptic time 
constants x^ and x,, similar to the VmD method (see discussion in Piwkowska et al., 2008). Interestingly, for 
some cases, the estimation works quite well (see indexed symbols in Fig.©. In the pool of injections, there are 
three cases that represent a cell in a high-conductance state, i.e. the mean inhibitory conductances are roughly 
three times greater than the excitatory ones, and the standard deviations obey a similar ratio. For these trials, the 
estimate comes close to the values used during the experiment. Indeed, we found that the relative error on a,- is 
roughly proportional to the ratio Oe/ot for ratios smaller than 1, and tends to saturate for larger ratios (Fig.O. In 
other words, the estimation has the lowest errors when inhibitory fluctuations dominate excitatory fluctuations. 
A recent estimate of conductance variances in cortical neurons of awake cats reported that a, is larger than Gg 
for the vast majority of cells analyzed (Rudolph et al., 2007). The same was also true for anesthetized states 
(Rudolph et al., 2005), suggesting that the VmT method should give acceptable errors in practical situations in 
vivo. 



Discussion 



In this paper, we have introduced and tested a new method to estimate conductance mean and variance from 
single V,,, traces. Traditional methods to estimate conductances require to perform cunent-clamp or voltage- 
clamp recordings at different voltages (Rudolph et al., 2004; Monier et al., 2008). If the membrane has intrinsic 
voltage-dependent conductances, the use of different voltage levels will necessarily cause important errors. 
This is particularly critical in vivo, where pharmacological manipulations to suppress currents are usually very 
limited. 

Several methods exists that can be applied to single V,„ traces; for example, one can calculate power spectral 
densities (PSD) of the V,„. The fitting of analytical expressions of the V„, PSD would allow in principle to 
estimate the conductance parameters, since they appear in these expressions. However, this analysis would 
require to know the effective membrane time constant, and therefore the total conductance. In such a case, 
since the average V,„ is known, the conductances can be estimated by simply solving the membrane equation at 
steady-state, so the PSD presents no useful means of estimating conductances. Nevertheless, fitting theoretical 
expressions to experimental PSDs can in principle yield estimates of other parameters, such as the decay time 
constants of synaptic conductances, although such estimates can be subject to important errors (see Piwkowska 
et al., 2008 for an assessment of this method). 

A second type of method applicable to single V„, traces is to estimate the spike-triggered average (STA) of 
synaptic conductances (Pospischil et al., 2007). The principle of this method is to search for the most probable 
conductance traces that account for the V,„ STA. These traces ai^e calculated by maximum likelihood following 
a discretization of the time axis. The performance of this method was tested using computational models and 
dynamic-clamp experiments (Pospischil et al., 2007; Piwkowska et al., 2008). It was found to be excellent, and 
the STA method was applied to estimate the optimal conductance patterns preceding spikes in vivo (Rudolph et 
al., 2007). However, this method requires a prior knowledge of the conductance parameters {geo, gm, Og, a,). In 
practice, it is thus necessary to apply the VmD method beforehand, which in turn requires recording at several 
V„, levels. 

We proposed here a method that merges the above concepts. The VmT method can estimate the same conduc- 
tance parameters as the VmD method (g^o, gjo, a^, Oi), but using a maximum likelihood approach similar to 
the STA estimates. The analysis can thus be realized from single V,„ levels, thereby limiting the bias caused 
by voltage-dependent cuiTcnts. The VmT method does not rely on static properties of the membrane voltage, 
like mean and standai^d deviation, in order to determine the shape of the conductance distributions. Instead, 
it exploits the dynamical information hidden in the V^ time course. Not only the step sizes At are impor- 
tant, but also at which voltage level they occur. The likelihood function is highly sensitive to the conductance 
fluctuations and, constraining the total conductance, also to the conductance mean values. 

A word of caution is needed here about the biological significance of the conductances estimated by the VmT 
method. Like any other method to extract conductances from V,„ activity, the conductance estimated here 
refers to the conductance visible at the site of the V,„ recording, in general the soma. This somatic conductance 
can be very different - and is in general smaller - than the conductance present at the site of the synaptic 
receptors, mostly the dendrites. This is especially critical for neurons with extended dendritic morphologies, 
such as pyramidal neurons, which experience considerable dendritic attenuation. However, the visible somatic 
conductances are close to the conductances participating to the interactions involved in generating the action 
potential, because the region for spike initiation (presumably the initial part of the axon; see review by Stuart 
et al., 1997) is electrotonically close to the soma. To yield estimates closer to the conductances present in 
dendrites, recordings must be done using drugs such as QX-3 14 and cesium, to reduce the leak conductance and 
thereby lead to more compact neurons, with more limited dendritic attenuation. The conductance contribution 
due to spikes (somatic or dendritic) should be avoided by using exclusively subthreshold activity. 



Similarly, like all other methods for the analysis of synaptic conductances, the VmT method rests on the as- 
sumption that the range of V,„ analyzed is linear, and does not contain voltage-dependent conductances. Nev- 
ertheless, it is possible to follow a similar maximum-likelihood procedure by explicitly including fast voltage- 
dependent conductances in the model. The VmT method also requires knowledge of the passive parameters 
of the recorded neuron (Gl, C.) which may be difficult to obtain in vivo (see Piwkowska et al., 2008 for a 
discussion of this point). 

A second word of caution concerns the significance of conductance variances. The conductance variances 
estimated using the VmT method are specific to the "point-conductance" model, in which the V,„ activity is 
assumed to result from the action of two stochastic synaptic conductances (Destexhe et al., 2001). In particular, 
one of the assumptions was that the excitatory and inhibitory conductances are Gaussian-distributed. Estimates 
of variances of course depend on this assumption. Including asymmetric conductance distributions is possible 
but would substantially complicate the model. The mathematical simplicity of this model has made it possible 
to design a series of methods to analyze the V„, activity (reviewed in Piwkowska et al., 2008). The VmT method 
is a natural extension of this approach. 

The tests of the VmT method on model data were encouraging, but also pointed to weaknesses of the method in 
a regime where the transmembrane current due to inhibitory conductances is small compared to the leak current. 
In this case, the inhibitory fluctuations are not reliably resolved. Furthermore, the presence of recording noise 
can disturb the estimation significantly. Though smoothing the voltage trace could in some cases re-establish a 
good result (Fig. O, more sophisticated procedures should be considered to eliminate this sensitivity to noise. 
For example, one could introduce to Eq. [T] an additive noise term, whose SD would become an additional 
parameter that has to be estimated. Thus, in principle, if this weak inhibition regime can be avoided, the VmT 
method should be applicable to real neurons. In any case, to be on the safe side, it is preferable to use the VmT 
method on V„, traces with the lowest possible level of instrumental noise. 

We also tested the VmT method by using dynamic-clamp injection of stochastic conductances. In this case, 
controlled conductances are injected while the V,„ activity is being recorded, allowing to compare the VmT 
estimates from the V„,, to the conductances actually injected. These tests confirmed the weakness of the method 
in the low-inhibition regime, but also revealed a good performance of the VmT method during HC states. Since 
in cortical neurons in vivo the membrane potential is usually well above the inhibitory reversal potential, and 
since inhibitory conductances tend to dominate over excitatory ones (Borg-Graham et al., 1998; Rudolph et al., 
2007; note that dendritic recordings may be dominated by excitatory inputs), we anticipate that this method 
should provide good conductance estimates from in vivo recordings. 
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Figure 1 : Test of the single-trace VmT method using a leaky integrator model. 
Each panel presents a scan in the {geo,gio)-plane. Color codes the relative deviation between model 
parameters and their estimates using the method (note the different scales for means/SDs). The white 
areas indicate regions where the mismatch was larger, >5% for the means and <-25% for SDs. A. 
Deviation in the mean of excitatory conductance (^eo)- B. Same as A, but for inhibition. C. Deviation 
in the SD of excitatory conductance. D. Same as C, but for inhibition. In general the method works 
well, except for a small band for the inhibitory SD. 
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Figure 2: The estimation error depends on the ratio of inhibitory and leak conductances. 
The relative deviation between the parameter a,- in the simulations and its re-estimated value is shown 
as a function of the ratio of the currents due to inhibitory and leak conductances. The estimation fails 
when the inhibitory component becomes too small. The same data as in Fig. [Hare plotted: different 
dots correspond to different pairs of excitatory and inhibitory conductances (geo, gjo), and the dots are 
colored according to the excitatory conductance (see scale). Details about the duration of injections 
etc. are given in the text. 
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Figure 3: Error of the VmT estimates following addition of white noise to the voltage trace. 
Gaussian white noise was added to the voltage trace of the model, and the VmT method was applied 
to the Yin trace obtained with noise, to yield estimates of conductances and variances. Left: relative 
error obtained in the estimation of geo and gio (upper panel), as well as Og and a,- (lower panels). 
Right: same estimation, but the V,„ was smoothed prior to the VmT estimate (Gaussian filter with SD 
of 1 data point). In both cases, the relative error is shown as a function of the white noise amplitude. 
Different curves correspond to different pairs (^eo^ ,?/o)- The errors on the estimates for both mean 
conductance and SD increase with the noise. The coloring of the curves as a function of the total 
conductance (see scale) shows that the largest errors generally occur for the low-conductance regimes. 
The error was greatly diminished by smoothing (right panels). 
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Figure 4: Dynamic-clamp test of the VmT method to extract conductances from guinea-pig visual 
cortical neurons in vitro. 

Fluctuating conductances of known parameters were injected in different neurons using dynamic- 
clamp, and the ¥,„ activity produced was analyzed using the single-trace VmT method. Each plot 
represents the different conductance parameters extracted from the V,„ activity: ^eo (A), g/o (B), Og 
(C) and O/ (D). The extracted parameter (Estimated) is compared to the value used in the conductance 
injection (Control). While in general the mean conductances are matched very well, the estimated SDs 
show a large spread around the target values. Nevertheless, during states dominated by inhibition (see 
indexed symbols), the estimation was acceptable. 
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Figure 5: Relative error on inhibitory variance is high only when excitatory fluctuations dominate. 
The relative mean-square error on a, is represented as a function of the Og/Oi ratio. The error is 
approximately proportional to the ratio of variances. The same data as in Fig. |4] were used. 



